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Abstract 



We derive a product rule for gauge invariant Weyl symbols which provides 
a generalization of the well-known Moyal formula to the case of non-vanishing 
electromagnetic fields. Applying our result to the guiding center problem we 
expand the guiding center Hamiltonian into an asymptotic power series with 
respect to both Planck's constant h and an adiabaticity parameter already 
present in the classical theory. This expansion is used to determine the influ- 
ence of quantum mechanical effects on guiding center motion. 

PACS: 03.65.Sq, 03.65.Ca, 41.75-i 



Typeset using REVT^X 



* email: mmueller @ eckar t . Ibl . gov 

^Present address: M. Miiller, Petrinistr. 7, 97080 Wiirzburg 

1 



1. INTRODUCTION 



In many physical applications charged particles are exposed to strong time-independent 
magnetic fields B{x) and additional electrostatic potentials 0(cc). Important examples are 
the magnetic confinement of plasmas, trapping of ions in accelerator facilities [Q as well 
as the Quantum Hall effect 0. Classically, the motion in such field configurations may be 
visualized as a fast rotation in the plane perpendicular to the magnetic field ("gyration"), 
with the center of the circular orbit moving slowly parallel to the magnetic field lines and 
drifting very slowly across both electric and magnetic field lines ("guiding center motion"). 
The underlying assumption that clearly distinguishable time scales of motion exist is known 
as "guiding center approximation" or, more general, "adiabatic approximation." 

To separate gyration and guiding center motion, different kinds of perturbative calcula- 
tions have been applied in classical mechanics 0-0]. Adiabatic invariants and equations of 
motion for the guiding center may be derived in a systematic way from Hamiltonian theory 
Pl^rSI . For a semiclassical description of guiding center motion the method invented by Lit- 



tlejohn [10,11 1 using non-canonical, but gauge invariant phase space coordinates turns out 



to be the best starting point. There, like in all classical investigations of the guiding center 
problem, a dimensionless expansion parameter e is introduced by replacing the electric charge 
q with q/e @,|T3- Physically, e represents the ratio of the gyroradius to the scale lengths 
of the external fields and is interpreted as an adiabatic parameter. Employing symplectic 
geometrical techniques, relations between the guiding center (phase space) coordinates and 
the particle's position and velocity are obtained which take the form of asymptotic power 
series in e. After writing down the Hamiltonian in terms of the guiding center coordinates, 
its dependence on the rapidly oscillating gyration angle is removed by means of averaging 
Lie transforms. The equations of motion resulting from the guiding center Hamiltonian 
confirm that the magnetic moment caused by the gyration is an adiabatic invariant. 

In low-temperature experiments the total energies of the particles are of order of the 
lowest Landau levels in the magnetic field B{x). Therefore quantum mechanical effects 
have to be taken into account when deriving equations of motion for the guiding center. 
So far this has been done only in the special case of a charged particle in the magnetic 
field outside of a rectilinear current filament |[T^|1^. To determine explicitly the quantum 
corrections to guiding center motion in arbitrary field configurations a method is needed 
which results in an expansion of the quantized guiding center Hamiltonian into a formal 
power series in both the (classical) parameter e and Planck's constant h. The quantum 
guiding center theory developed by Maraner in two inspiring papers p!7| , p!8[] uses only the 



magnetic length Ib = yhc/ {q\B\) as expansion parameter. The power series expansion of the 
guiding center Hamiltonian operator with respect to Ib does not distinguish between terms 
of adiabatic origin already present in classical mechanics and quantum corrections caused 
by the non-commutativity of the operator algebra. For experimental purposes, however, it 
is very important to know whether the classical picture is valid even at low temperatures or 
whether quantum effects dominate guiding center motion. 

A first step to answer this question lies in the observation that there is great formal 
resemblance between guiding center motion and adiabatic motion of neutral spinning parti- 
cles in an inhomogeneous magnetic field |]19[. The latter has recently been studied in more 
detail because it represents a standard example for the occurrence of "geometrical" forces 



2 



in dynamical systems pO|-p^. A semiclassical investigation of this motion p4| involves a 



multicomponent version of the Weyl calculus p5|-p7| . It has the appealing feature that two 
different expansion parameters are used: one, e^, connected with adiabaticity (i.e. the as- 
sumption that the magnetic field does not change appreciably during a precession period) 
and another, e^, proportional to controlling the validity of the semiclassical approxima- 
tion. In the diagonalized Hamiltonian, which describes orbital motion, the potential terms 
are expanded with respect to both ea and e^. To achieve the same goal for the guiding center 
Hamiltonian, the Wigner-Weyl formalism |2^j2^ therefore seems to provide the appropriate 
tools. 

In general the Weyl transform of a quantum mechanical operator is a uniquely determined 
phase space function which may be defined as follows pOH3^: Starting from the fundamental 
operators x and p, a particular, continuously indexed basis A(a;,p) of the operator space 
is constructed. (Here, as in the following, the hat denotes an operator.) The representation 
of an arbitrary operator A as a linear combination of the operators A(a;,p) involves c- 
number coefficients which are labeled by the continuous variables x and p. They constitute 
a function Aw{x, p) on phase space which is denoted as the Weyl symbol of the operator A. 
The relation between the symbol Cw{,x,p) associated with an operator product C = AB 
and the symbols Aw[x,p) and Bw[x,p) of its factors is given by a nontrivial composition 
rule known as Moyal formula p3| . 

If a magnetic field B{x) = S/xA{x) is present, however, the Weyl correspondence should 
be re-defined, because the gauge dependence of the canonical momentum p causes the basic 
operators A{x,p) to be gauge dependent as well. As shown in |Q this leads to the un- 
desirable consequence that the Weyl symbol of a gauge invariant operator becomes gauge 
dependent and that, vice versa, the operator corresponding to a gauge invariant phase space 
function is itself in general not gauge independent. The most natural way to include the 
principle of gauge invariance into the Weyl formalism is to replace the (gauge dependent) 
canonical momentum p appearing in the definition of A{x,p) by the (gauge invariant) ki- 
netic momentum k = mv = p — {q/c)A{x) ||3^-p^. The coefficient function of an operator 
A (not to be confused with the vector potential) with respect to the new set of basic opera- 
tors A{x, k) will be denoted in the following as the gauge invariant Weyl symbol Aw[x, k) 
of A. One can show that Aw[x, k) is a gauge invariant phase space function if and only if 
A is gauge invariant (for more details cf. and section 2). Obviously, the product rule 
for gauge invariant Weyl symbols will be different from the usual Moyal formula. 

After comparing ordinary and gauge invariant Weyl calculus in a little more detail, we 
will explicitly derive the gauge invariant generalization of Moyal's formula in the next section 
and discuss its most important properties. In section 3, the gauge invariant Weyl formalism 
will be applied to separate the different time scales occurring in the motion of charged 
particles in external electromagnetic fields within a semiclassical framework. As a result 
we will expand both guiding center coordinates and the guiding center Hamiltonian into 
asymptotic power series with respect to the adiabatic parameter e and Planck's constant 
h. Section 4 contains a summary of our results and conclusions concerning the infiuence 
of quantum effects on guiding center motion which can be derived from our expansion of 
the guiding center Hamiltonian. Finally we compare our results to the quantum mechanical 
calculations of Maraner. In the appendix, the classical guiding center theory for the motion 
of a charged particle in a magnetic field of constant direction and an additional electrostatic 
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field is summarized. 



2. PRODUCT RULE FOR GAUGE INVARIANT WEYL SYMBOLS 



In order to set the stage for our computations, let us briefly review some basic features 
of the ordinary Weyl transform valid in a six-dimensional flat phase space in the absence of 



magnetic fields ||30 H32| , |37| , p8| . Starting with the set of generating Heisenberg operators 

T{u,v) = exp[i{u-p + v-x)] , (1) 



we introduce a basis 



— ^ / d'^u d^v exp[i{u-p + v-x)] T{~u, —v) 
2^ 



I j d^ud\eyi:{i{i\u-{p — p) + v-{x — x)\} 



(2) 



of the operator space which is labeled by the continuous classical variables x and p. If an 
operator A is written as a linear combination of the A(a;,p), 



A 



d^x d^ 



A{x,p)Aw{x,p) , 



(3) 



(4) 



is called the Weyl symbol associated with A. Note that equation (|) is a direct consequence 
of definition fRI) and 



the uniquely determined coefficient function 

Aw{x,p) = fd^x'{x'\AA{x,p)\x') =Tr[iA(x,p) 



Tr 



A{x, p) A{x\ p') = {2nhf6{p - p') 5{x - x') . 



(5) 



Here, as in the following, we leave aside questions of convergence and the mathematical 
problem of characterizing the class of operators for which expansions like @ exist. 
From equations (|T|)-(^) and the duplication formula 



f{u,v)f{u\v') = exp 



i — (u-v' — u'-v] 

2 ^ ^ 



(6) 



one can immediately determine the relation between the Weyl symbol Cw{x,p) of a product 
operator C = AB and the symbols of its factors. The result is the well-known Moyal formula 



Cvy(a;,p) = exp 



ihJ^/d d _ d d ' 
2 ^^\dbi dzi dtti dyi^ 



Awiz,a)Bwiy,b) 



y^z^x 
a=h=p 



(7) 



where the subscript "i" characterizes the Cartesian coordinates of a vector, i.e. 1,2,3 stands 
for x, z respectively, and the auxiliary vectors a, 6, t/, z specify which of the factors 
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Aw{x, p), Bw{x, p) has to be differentiated with respect to x and p. Note that the operator 
in the exponential is just the ordinary {x,p) Poisson bracket operator, so that expanding 
the right hand side of (|^ with respect to h yields 

Cwix,p) = Awix,p)Bwix,p) + '-^{Aw,Bw} + 0{h^) . (8) 

Equation (|^) may also be interpreted as defining a bilinear, associative and non- commutative 
product on the space of symbols, 

Cw{x, p) = Aw{x, p) * Bw{x, p) , (9) 

denoted as the star product or Weyl product. 

Suppose we are given a phase space function of the form 

f{x,p)=xTp], (10) 

with m,n & INq (= non-negative integers). Evaluating equation we obtain the corre- 
sponding operator |^ 

1 ™ /'m\ 

/(^>p) = ^E( J^iPi^r'- (11) 

It can be constructed in the following way: First take Xi m times, pj n times, put them 
in all possible permutations with equal weights and divide by the number of terms. The 
result is called the totally symmetrized or Weyl ordered product, written as Symm^f^p"^. 
Finally apply the commutation relation [xj,^^] = ih6ij to bring the p/s together at various 
positions of the product with no terms proportional to h remaining. Due to the linearity of 
the Weyl transform, equations (^) and (^) generalize to any analytic function on phase 
space. 

So far we have used position x and canonical momentum p as basic variables. If a mag- 
netic field B{x) = VxA{x) is switched on, the canonical momentum is no longer completely 
physical because of its gauge dependence. However, the operator of the kinetic momentum, 

k = p A{x) = mv , (12) 

is gauge invariant because its expectation value is not effected by a gauge transformation 
40| . In contrast to p the Cartesian components of k do not commute with one another. 



q 

ih-eijkBk{x), (13) 



whereas their commutation relations with x parallel those for p and x, [xi, kj] = ih6ij. 

If we replace the canonical momentum p in (|l|) with the kinetic momentum k, the 
operators 

T{u,v) = exp i(^u-k + v-ihj (14) 
become gauge invariant and hence the basic operators 
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2tt 



d^u d^v exp[i{u- k + v-x)] T{—u, —v) 



(15) 



are gauge invariant as well. The Weyl symbol Aw{x,k) of an operator A is now defined 
with respect to the new basis in the same way as in the field-free case, 



Awix,k) = Tr AA{x,k) 



or, equivalently. 



A 



d^x d% 



A(a;, k) Aw{x, k) 



(16) 



(17) 



From equation ([T6|) and the properties of gauge invariant operators (cf. |0) it is obvious 
that the Weyl symbol of a gauge invariant operator does not change its value under gauge 
transformations. Thus, viewed as a phase space function, the symbol is also gauge invariant. 
According to equation (|1^) the opposite is also true: The operator corresponding to a gauge 
invariant phase space function is itself gauge invariant, i.e. its mean value does not change 
under gauge transformations. 

Writing Aw{x, k) Fourier integral. 



Aw{x, k) = d^ud^v exp[i{u-k + v-x)] A{u, v) , 



(18) 



one can show by inserting definition ([T5| ) into (|T6D that the Fourier transform A{u, v) may 
also be obtained from 



A{u,v) 



2tt 



Tr 



Af{-u,-v] 



(19) 



According to (|T3|), (0), and (|TB|) the operator A can similarly be expressed in terms of 
A(u, v) via 



A= Id^u d^vA{u, v) T{u, v) 



(20) 



Equations (|14|), (|T8|) and (^) are the gauge invariant generalization of Weyl's original cor- 
respondence rule for phase space functions and quantum mechanical operators. 

To derive a product rule for gauge invariant Weyl symbols we will have to evaluate matrix 
elements of the form (a;'|yl|a;). For this purpose it is of advantage to express the operator 
exp(iii-fc) occurring in T{u,v) by the translation operator exp{iu-p). The latter acts on 
an operator function f{x) and a position eigenstate la;) in the following way 



exp{iu-p)f{x) = f{x + hu) exp(i u-p) 
exp(i u p) \x) = \x — hu). 



(21) 
(22) 



As proven in pB|, exp(?tt-fc) and exp(i'u-p) are related by 
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exp(i u-k) = expU up) exp| —i-u- A(x — %tu) dr 

I C JO 

= expi —i-u- A{x + fiTu) dr i exp{iu-p) . 



(23) 



Using (|2l|), (|23|) , and the Baker- Campbell-Hausdorff type formula 



T{u, v) = exp(— i u-v/2) exp(i u-k) exp{iv-x) 
= exp(i u-v/2) exp{iv-x) exp(i u-k) , 



(24) 



a straightforward calculation shows that the product of two gauge invariant T operators 
may be cast into the form 



t(u, v) f(u', v') = exp<{ i-u-[ \A(x + firiu + u')) - Aix + Titu) 

c Jo L 

r-l 



dr \ X 



exp < 2 - w ■ 

' c Jo 



exp 



z — n-t; — n -f 

2 ^ ' 



A{x + /ir(w + u')) - A{x + hu + hru'] 
f{u + u\v + v'), 



dr > X 



(25) 



which reduces to the ordinary duplication formula (|^) if the vector potential A{x) vanishes. 
Another consequence of equations ( p3D and ( p^ is that the trace of T(w, v) T{—u', ~v') is 
given by 



Tr 



T(ii, v) T{—u', ~v') 



2tx 



b{u — u') 5{v — v') 



(26) 



After these preliminary remarks we are ready to determine the Weyl symbol Cw{x,k) 
of the operator product C = AB in terms of Aw{x, k) and Bw{x, k). From equations 
and (|25|) we find that 



C = AB= I d^u d^u' d^v'A{u, v) B{u\ v') T{u, v) T{u', v' 
d^u d^v d^U d^VA{u, v) B{U - u,V - v)x 



exp 



i- (u-V -U-v) 

2 ^ ^ 



F{u,U] x)f{U,V) 



(27) 



where we introduced new integration variables U=u + u', V = v + v' and the function 
F{ui, U2] x) = exp/z - Ui- [ [ A{x + h{l — t)ui + hTU2) — A{x + hrui)] dr — 

I C JO 

i-U2- [ A{x + h{l — t)ui + hTU2) — A{x + fiTU2)] dr \ 

C JO J 



exp 



i- f{ui,U2;x) 



(28) 



which is equal to unity in the field-free case A{x) = 0. Note that in (^) F depends on the 
position operator x and hence is itself an operator. 



7 



Inserting p7| ) into (|T9D and making use of (^) and the completeness of the position 
eigenstates \x) yields the Fourier transform of Cw{x, k), 

C{u,v)= l^^j J d^xd^u'd^v'd^VA{u',v')B{u-u',V -v')exp[i{V -v)-x]x 

expj i ^ [{u + u')-V -u-{v + v')] I F{u', u; x) , (29) 

so that as an intermediate result the Weyl symbol of C reads 

Cw{x, k)= [ d^u' d^v' d^u" d^v"A{u', v') B{u'\ v") exp{ t [{v' + v") ■x + {u' + u") ■k]}x 



exp 



h 

% - iii! -v" - v!' -v') 
2 ^ ' 



(30) 



where u" = u — u', v" = V — v' . As a next step we want to express the right hand side of 
(pop in terms of Ay/{x, k), B]y[x, k) and their derivatives with respect to x and k. A helpful 
observation is that by setting A{x) = the above equation reduces to the one appearing in 
the derivation of the ordinary product rule ||S^. There, the factor exp[ih {u' -v" — u" ■v')/2] 
is expanded into a power series with respect to %. The variables Ui,u'- ,vl,v" occur- 
ring in each term of this series are generated from exp{i [{v' + v")-x + {u' + u")-k]} = 
exp[i{v' -x + u' ■k)]exp[i{v" -x + u" -k)] by differentiation processes. This is achieved by 
replacing the variables x and k with z and a in the first exponential factor and with y and 
b in the second one and then applying the operator ih{d/dz-d/db — d/dy-d/da)/2 and 
appropriate powers of it to the product. The resulting total differential operator has the 
form exp[i^ [d/dz-d/db ~ d/dy-d/da)/2]. If it is taken outside of the integral, the latter 
may be evaluated and one finally gets Moyal's formula (^. 

To employ a comparable algorithm in the case of non-zero vector potential we have to 
extract u' and u" from the integrals in F{u', u' + u"\ x — h{u' + u")/2). For this purpose 
the vector potentials occurring in the exponent of (^) are expanded into Taylor series around 
the position x. After some additional algebraic manipulations we obtain 



f{u\ u'+ u"- x-h{u' + u")/2) = E ^ E 



n=0 



..in d- 



- < I n [ (1/2 - r) - 1/2 <J dr + u'l H [ V2 + (1/2 - r) 

n I 



dr 



(31) 



Noting that for /c G INc 



l/2)'=rfr 



1\ fc+l _ 

2J k + l 



if = 2m + 1, me INq, 



i) ^ > , if = 2m, m G INq, 



(32) 
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we arrive at 



f{u', u + u"- x-n{u' + u") /2) = 



n=l 



71+1 



{n + iy 



Ir 



X 



r,j,l. 



u 



fc=i 



, k , 



(i-(-i)') (n+i)-(i-(-ir^)fc 



dxi^ . . . Oxi^_^ 

(33) 



The most important result of the preceding calculation is that the expansion of / includes 
only derivatives of the magnetic field B{x). Physically this was to be expected because the 
product of two gauge invariant symbols is itself gauge invariant. Therefore the integrand on 
the right hand side of ( pOf ) must not depend on the chosen gauge. As all other factors satisfy 
this condition, the function F{u', u' + u"; x — h {u' + u")/2) has to be gauge invariant as 
well. This is certainly true if it is a functional of the magnetic field. 

Now we continue just like in the field-free case. Writing 

exp{ i[{v' + v")-x + {u' + u")-k]} = exp [i{v'-z + u'-a)] exp[i{v"-y + u"-b)]\y=z=x (34) 

we may generate the variables m^, u", f^, v" occurring in an analytic function which is multi- 
plied to the right by this exponential by differentiating the latter with respect to the auxiliary 
variables ai,bi,yi, Zi. This is formally equivalent to substituting 



d 



■ 8 , 



d 



d_ 



(35) 



in the analytic function. Hence, if we introduce gauge invariant operators 



~ 2 dai dyi dk dzi ' 



(36) 



n+1 



E 



^jlr 



^u. 2^ 



{n + 1)2 n\ f-'i^ dxi^ . . . <9xi„_, daj dbi k 

n,...,i„_i=i 



(n + 1) - 1 



-1 ' 



n+l 



d 



d d 



X 



d 



(37) 



G IN, and leave aside questions of convergence, we may write 

exp [ i (n' ■ v" - u" ■ v') /2] exp{ i[{v' + v")-x + {u' + u")-k]} 

exp[— ey.-p[i{v' -z + u'-a)] exp[?(v"-j/ + u" •b)]\y^z=x , 

a=b^k 



(38) 



and 



F(n', u' + u"- x-n{u' + u") /2) exp{ i [{v' + v") • a; + (n' + u") ■ k] } 



exp 



^ E ^'^^n 



n=l 



exp[i{v'-z + u'-a)]exp[i{v"-y + u" ■b)]\y=z=x ■ 

a=b=k 



(39) 
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All operators £ and n ^ IN, commute with one another because they contain 
ai,bi,yi,Zi only as differentiating variables (the magnetic field B{x) and its derivatives 
occurring in depend on the position x and are hence not effected by a differentiation 
with respect to these variables). Therefore, the integrand in (^) can be generated by the 
action of 



V = exp 



-ihc-i{q/c)Y,n''Cn 



n=l 



(40) 



on exp [i{v'-z + u' -a)] exp[i(i7"-?/ + u" ■h)\\y^z=x ■ Taking the total differential operator 

a=b=k 

outside of the integral (|30D we finally obtain 



Cw{x^ k) = exp 



n=l 



Aw{z,a)Bw{y,b) 



y^z=x 
a=b=k 



= [Aw * Bw]{x,k) 



(41) 



which is the generalization of Moyal's formula to gauge invariant Weyl symbols. As in the 
case of the ordinary Weyl transform one can show that the star product defined by ( ^I]) is 
bilinear and associative. 

Before turning to the semiclassical analysis of guiding center motion, let us investigate 
equation (|4l|) in more detail. Expanding the exponential operator V into a power series 
with respect to h yields 



V =l-ih{C + ^Ci] -h' 



(42) 



The second term on the right hand side of ( ^2] ) is equal to ih/2 times the {x,p) Poisson 
bracket operator. This is seen most easily by expressing the Poisson bracket of two arbitrary 
phase space functions in terms of x- and fc-derivatives. 



jr ._^df dg df dg q df dg 



~^ dxi dki dki dxi 



dkj dki 



(43) 



A comparison of (^ and (^) shows that the first order terms of both expansions coincide. 
However, the higher order terms in (P) turn out to be gauge dependent and hence differ 
from those in (|4^) . 

Using their definitions (pG]), (p7D, one can derive the following symmetry properties of 
the operators C, Cn, 



/:^{Bw,Aw) = (-1)"'"C(Ah.,5^), meIN 



(44) 
(45) 



which cause the star product to be non-commutative. The difference 
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Aw*Bw -Bw*Aw = [Aw, Bw]m (46) 
is called the Moyal bracket of Aw and Bw- According to (|4^ ) it may be expanded into 



[Aw J Bw]m = i^wBw) 



c \ c J c 



+ 0{h') AwB 



(47) 



where the leading order term is just % h {Aw, Bw}- In contrast to the ordinary Weyl calculus 
the gauge invariant Moyal bracket operator Ai cannot be written in closed form. 

Finally one can prove from (^T]) by induction that the Weyl symbol of the operator 



is given by 



A = Symml^ki.ki^ . . . , ij E {1, 2, 3}, 1 <j <n, 



Aw = fc3^ m + 712 + ng = n , 



(48) 



(49) 



if fcj, 1 < i < 3, appears rij times in the operator product kij^i^ ■ ■ - kt^. The relation above is a 
direct consequence of the non-commutativity of the Cartesian components of k. In general, 
if f{x) is an analytic function of the operator related to 



Aw = fix) A;^ k^ k 



713 
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reads 



1=0 



A = Symm^^^ J (^] ■ ■ ■ K-i /(*) K-i+iK-i+i ■■■kin 



(50) 



(51) 



where n = ni + ?t,2 + ?t-3 and Symm|^ ^ denotes symmetrization with respect to the operators 
kj - . 



3. SEMICLASSICAL DESCRIPTION OF GUIDING CENTER MOTION 

We will now apply the gauge invariant Weyl formalism to describe the motion of a 
charged particle of mass m in a strong time-independent magnetic field B{x) = VxA{x) 
and an additional electrostatic field E{x) = — V0(cc) semiclassically. The influence of the 
latter has not been taken into account in quantum mechanical calculations ||l7| , p!8[| so far. 



Assuming that the guiding center approximation is valid, we are in particular interested in 
the lowest order quantum mechanical correction to the guiding center Hamiltonian. 

To incorporate the guiding center approximation into our theory, we follow the classical 
calculations and introduce an adiabatic parameter e by replacing the electric charge q of the 
particle by q/e [|TU],|TT|, 



g - f ■ (52) 
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Physical results are recovered at the end of our calculation by setting e = 1. In guiding center 
approximation we are speaking of the order of an expression in terms of its behavior as e— >0. 
The physical meaning and mathematical details of this limit are discussed in greater detail 
in PJI^. We adopt the convention that the particle variables x and v as well as the fields A 
and B are held constant in this limiting process, i.e. are independent of e. Since the guiding 
center approximation breaks down when the component E\\ of the electric field parallel to 
B is of the same magnitude as \B\, we take E\\ = 0{e) []TT|. 

For reasons of notational convenience we will further on suppress the constants g, m, c, 
which is equivalent to the following scaling of physical quantities, 



1 



X 



X, p 



mp, V 



m 



"^'a^b-^^b. 



E 



m 



E. 



(53) 



The corresponding backward transformations restore the correct physical units in our results. 
Note that with respect to this scaling particle velocity and kinetic momentum are equal. 
Due to the foregoing conventions, the operator V introduced in (^) takes the form 



V = exp 



n=l 



(54) 



where £, £„ denote the scaled versions of Taking into account the definitions (p6D, 

(^) of £, Cn the first three terms of the power series expansion of V with respect to h read 



E 



1 ' 



/I I o E 



■[ \daidyi dhdzi 



1 ^ 



dajdbi 



+ 



4 I 2 f-^-^ \daidajdyidyj daidbjdyidzj dbidbjdzidz 



1 3 



J — I y J-JO-JJ 1,- - ct- -J ~ - c~ -J - -0- -J / 

9^ \ ^ 1 dBr ( 



Q3 



dttidajObidyi dajdbidbidzi J 3 dxi \dajdaidbi dajdbidb 



1 1 3 

+ 2^ '^"-'""^ da.dakdbidbr, 

fc.m.s — 1 



Oh 



(55) 



Since each of the operators has attached to it, a simple reasoning shows that the 
term proportional to ^" in the expansion of V includes terms of all orders in from to 
n. Therefore we may formally write 



(56) 



n=0 



n=0 m= 



The operators Vn,m contain derivatives of order n— |m| with respect to the position variables 
y,z and derivatives of order ?T,+ |m| with respect to the kinetic momentum (=velocity) 
variables a, b. As a consequence of the symmetry properties of C, Cn no terms of even 
power in h occur in the expansion of the (scaled) Moyal bracket operator Ai, 
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M = 2j2h'''+'V2n+l = 2Y: E ^'"+'e'"P2n+l,n^. (57) 
n=0 n=0 m=-(2n+l) 

Suppose now we are given symbols of the form 

oo oo oo oo 

Aw = = E Bw = E^"^n = E , (58) 

n=0 m,n=0 n=0 •m,n=0 



with the coefficients being analytical functions of x and k. From ( p6D and (|58|) we obtain 
the following expansion for the star product of Aw and Bw, 

oo n oo m+i n—i I 

Aw*Bw = Y.Y. E E E E ^r,s5n-.-r,i-s • (59) 

n=0 i=0 m=—i l=m r=0 s=0 

A similar calculation yields for the Moyal bracket of Aw and Bw 



[Aw, Bw]m — E ^'^ 

n=0 



E [l-{-ir]V^[ E AiB^_^_i 

.m=0 \ /=0 > 



oo 71 oo m+i n—i I 

= T.T. H ^"e™ [l - (-1)* ] V.,m-l E E Ar,sBn^^-r,l-s ■ (60) 

n=0 i=0 m=— i Z=?Tt r=0 s=0 

Let us now determine the symbols of both guiding center coordinates and the guiding 
center Hamiltonian. The basic features of our method become most transparent if the 
direction of the magnetic field is constant with the electric field being perpendicular to 
it. In addition, this case is notationally easier to handle than the more general one in 
which the directions of both fields are varying arbitrarily. Therefore, we will consider in the 
following a charged particle in a magnetic field B = B{x, y)ez and electrostatic potential 
(f){x, y) neglecting its motion parallel to B. This kind of planar motion in a strong two- 
dimensional magnetic field (i.e. of constant direction) is intensively studied in the context 
of the Quantum Hall effect 0. There, the electric field is weak compared to B, which in 
our scaling is equivalent to assuming that E is of order e, E = 0{e). Classically, this means 
that the ExB drift is of the same order of magnitude as the VB drift 0. Using scaled 
velocity operators Vi = pi — e~^Ai = ki, i = 1,2, the Hamiltonian in such a field configuration 
reads 



H 



l{vl + vl)+4^{x,y). (61) 



Its Weyl symbol is obtained by replacing operators with their corresponding phase space 
functions, taking into account that the symbol of vf = ViVi is equal to Vi*Vi, 



Hw = ^{v^*v.^ + Vy* Vy) + (f){x, y) . (62) 



In the special case of a homogeneous magnetic field B = Bez, it is well known that the 
operators 

V, = B-'/^v,, Vy = B-'/\ (63) 
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are canonically conjugate and the Hamiltonian has the form of a one-dimensional harmonic 
oscillator. Physically, and Vy describe the gyration around the magnetic field lines. To 
get a complete set of conjugate operators including and Vy, one has to replace the particle 
coordinates with the operators 

X = x + eB-^Vy , Y = y- eB'^v^ (64) 

of the guiding center position. In equations (|63| ) and (|6^ ) questions of ordering need not be 
taken into consideration because 5 is a real valued constant. The non- vanishing commuta- 
tors of X, Y, Vx, Vy are 

[Vx,Vy]=z-^, [X,Y]=^^. (65) 

The Weyl symbols of these operators are obtained by replacing x,y,Vx,Vy in (|63[) and (|64D 
with the corresponding phase space functions. We denote them as guiding center symbols 
and - leaving away the subscript "W" - write X,Y,Vx,Vy or X,V for them. Their Moyal 
brackets resemble the commutators (|65D of the related operators. 

Generalizing the results for the homogeneous field to the case of an arbitrary two- 
dimensional magnetic field we are looking for a set of symbols {X, Y, Vx, Vy) whose non- 
vanishing Moyal brackets are given by 

[K.,v;]m = ^^, [^'^W = ^;^(^- (66) 

Their different orders with respect to e indicate the different time scales of motion. They 
are separated because the symbols X, Y of the guiding center position commute with those 
of the gyration velocity Vx, Vy. The latter are again conjugate to one another. 

Concerning the Moyal bracket of the guiding center position components X and Y two 
remarks are necessary: First, the symbol B{X, Y) specifies the strength of the magnetic field 
at the position of the guiding center. The corresponding operator is uniquely determined if 
X and Y are expressed in terms of the particle coordinates x, y, Vx, Vy and the correspondence 
rule (13) 5 ( |51| ) for arbitrary Weyl symbols and their operators is applied. Second, one could 
think of replacing X, Y with Euler potentials x^{X, Y),x'^{X, Y) thus obtaining a set of 
conjugate variables to describe guiding center motion. However, Euler potentials are non- 
physical in the same sense as the vector potential A is. Moreover, in a three-dimensional 
magnetic field we get four non-canonical guiding center coordinates instead of X and Y 



rT| . To transform them into two pairs of canonically conjugate variables one has to find 
functions which are less familiar than Euler potentials and much more difficult to construct. 
Therefore we will keep using non-canonical coordinates X, Y to specify the position of the 
guiding center. 

From classical guiding center theory ||10[ it is well known that {X, Y, Vx, Vy) can be chosen 
in such a way that J = Vx + Vy is a constant of motion which may be interpreted as the 
generalized magnetic moment of gyration. As a direct consequence of the relations ( |66D 
the guiding center Hamiltonian must therefore depend on Vx,Vy only by means of J and 
its powers. Accordingly, to find an appropriate set of guiding center symbols we have to 
proceed as follows: First we determine symbols X, Y, Vx, Vy satisfying (|66D. Next we express 
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the particle phase space coordinates {x, y, Vx, Vy) in terms of them and insert our result into 



the Hamiltonian (32). If the latter contains \4, Vy in other combinations than J we have to 
transform to a new set of averaged guiding center symbols {X, V) satisfying the same Moyal 
bracket relations, but Hw depending on the gyration velocities only via * Vx + Vy * Vy. 
This symbol transformation is an analog of the near-identity Lie transform carried out in 
the classical calculation [|TO| ,p]T[]. 

To begin with, let us analyze the Moyal bracket relations ( |5BD in more detail. Assuming 
that the guiding center symbols can be expanded into power series with respect to h and e as 
specified in (^), the Moyal brackets take the form (^) with almost all coefficients vanishing. 
Only those of he^^ and he are different from zero if A = V^, B = Vy and A = X, B = Y , 
respectively. As stated earlier, the ^-term of the Moyal bracket [Avi/, -Bvy]Af is proportional 
to the Poisson bracket of Aq and Bo. Therefore, the zero order terms (with respect to h) of 
the guiding center symbols satisfy the Poisson bracket relations of classical guiding center 
theory. Hence we identify them with the classical guiding center coordinates, in agreement 



with the fact that in the limit /i— >0 Weyl symbols become classical functions |3^. As we 
will refer to them frequently in the remainder of this section, the results of the classical 
guiding center theory in a two-dimensional magnetic field (using Cartesian coordinates) are 
briefiy summarized in the appendix. 

Concerning higher order terms of the /i-expansion of the guiding center symbols one can 
show that 

A2n+1 = , B2n+1 = , 71 G INq , (67) 

which means that only even powers in h occur. In addition, one can prove that the coefficients 
A2n,i are zero for < i < 2?7, — 1. For i>2n they turn out to be homogeneous polynomials 
of degree i — 2n + l for the components of V and of degree i — 2n for the components of X. 
Hence, the expansions of the guiding center symbols with respect to h and e take the form 



oo 

V= Y: h'-e'-+- Y: V,t'S+^{x, y) , (68) 

m,n=0 fei+fc2 

= m+l 



oo 



X,= ^'"e''^^"^ E <'vl^X}tZlrn{x.y). (69) 

= m 

where i = 1, 2 denotes the Cartesian components of V and X. The coefficient functions 
^2^n,'2n+m(^) y); -^/2n,'2n+m(3^? y) O)) (^1) fuuctiouals of the clcctric and magnetic 
fields and thus gauge invariant. 

To express the Hamiltonian (^) in terms of the guiding center symbols we have to find 
the corresponding backward transformations. Formally they are given by 

oo 

v,= Y ^'"^'""-"^ E V^^V^^v^t&^iX^y). (70) 

m,7i=0 fci+fe2 



oo 



Xi 



Y ^'"^'"^"^ E V^^V^'^x}tfJ^^{X,Y). (71) 

m,n=0 fci+fe2 



Again, the coefficient functions Vi2n]2n+m{X^y)^^i2n,2n+rn{X-,y) ^'^^ gauge independent. 
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The products of symbols in (|68|)-([7TD are defined point-wise. To derive corresponding 
relations between the particle operators x, i) and the guiding center operators X, V from 
(|68|)-([7T|), we have to introduce the star product on the right hand side of these equations. 
In (|68|) and (|^) this is simply done by making use of the correspondence rule 
In the case of the backward transformations (|70D, ([7I| ) we have to "translate" the classical 
result (|A9|) - (|A12|) given in the appendix into operator "language" to obtain the lowest order 



term with respect to h. It will turn out that this is enough to determine the lowest order 
quantum mechanical correction to the guiding center Hamiltonian. 

To this end we first conclude from equation ( |5BD that for any two symbols A^, Bw 
their point-wise product and star product differ by Y^'^=i^^'T^n{.A.w , Bw) . Since the Moyal 
brackets of the components (X, Y) of the guiding center position and the gyration velocity 
{Vxi Vy) vanish, we have for two arbitrary functions /(X, Y) and giVx-, Vy) 



fix, Y) * giVx, Vy) = giVx, Vy) * /(X, Y) (72) 



and 



P2n+i [fix, Y),giVx, Vy)] = 0, ne INq . (73) 

As an example, 

B'/\X, Y) * K. = K. * B'/^iX, Y) = B^/\X, + h^V^^B^'^ K) + O [h') , (74) 



with V2\B^/'^ ,Vxj being of order e, so that the difference between B^/'^V^ and B^/'^*Vx is 

of order h'^e. The fact that the symbols x,y,Vx,Vy represent self-adjoint operators leads 
directly to the substitutions 



V.Vy^liVx*Vy + Vy*V,), V^Vy ^ V, * Vy * V, , V^, ^ Vy * V, * Vy (75) 
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in equations (|A^)- (|XT2D , because the symbols on the right hand side of ( |75D correspond to 



self-adjoint operators. Further computations show that the replacement of the point-wise 
product with the star product in equations (|A9|)-( PT2D leads to corrections which are at 



least of order h^e or he^. Therefore, up to terms of order the relations between the particle 
operators (symbols) and the guiding center operators (symbols) formally coincide with the 
classical result ([X^)- (|XT2|) , if the point- wise product is replaced by the star product and 
the substitution rules ( [75| ) for products of gyration velocities are taken into account. For 
this reason we refrain from writing them down explicitly and refer the interested reader to 
the appendix. 

Next we insert the results for Vx,Vy into the Hamiltonian ( |62D and expand the potential 
0(x, y) into a Taylor series around the guiding center position (X, Y) replacing again point- 
wise products with star products. Naively one would expect from our previous results that 
the symbol Hamiltonian is of the form 

Hw = H,i + Oie%en^), (76) 

with Hci being formally equal (in the sense described above) to the classical Hamiltonian 
function ( [A13| ). However, as pointed out earlier, Hw should depend on the gyration velocity 
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components Vx, Vy only by means of the magnetic moment of gyration, J = Vx*Vx+Vy*Vy, and 
its powers. When evaluating the products Vx*Vx and Vy*Vy we have to change the ordering of 
the symbols Vx, Vy accordingly. Since their Moyal bracket is of order he~^, additional terms 
compared to the classical Hamiltonian occur. A straightforward calculation shows that the 
lowest order correction originates from the e^-term in the classical Hamiltonian function. 
It is of order h^. Leaving away from now on the multiplication symbol "*" for reasons of 
notational simplicity, the Weyl symbol of the guiding center Hamiltonian finally turns out 
to be 



Hw = lB{v,' + V;)+<l)iX,Y) 

(-Wl + BB^x. - ^Bl + BB^yy) ( + V; 
+ 4 {2>ExB^x - BEx,x + SEyB^y - BEy^y) {y^ + If) - 8 (eI + 



2 
+ 
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+ ^ {-Bl + BB^x. - Bl + BB^yy) + 0{eh\ e% e') . (77) 

Here, squares of symbols are star products of equal factors, electric and magnetic fields 
have to be evaluated at the guiding center position and the comma in a subscript denotes 
differentiation with respect to the following coordinate(s). Using two-dimensional vector 
notation equation ( |77D may be written in a more compact form, 

1 

Hw = -BJ + (f){X,Y) + —— \(bAB -3\VB\'^)j^ + A{3E-VB ~ BV-E) J -8\E 



+ ^ {BAB - I + 0{eh\ e% e') , (78) 

where scalar products denote summation over Weyl products of vector components. Due 
to the use of the gauge invariant Weyl calculus, the expansion of Hw involves only gauge 
invariant quantities. Thus, the higher order terms of the guiding center Hamiltonian will be 
gauge independent as well. 

Up to second order in h and e the Hamiltonian Hw in ([77|) depends on the gyration 
velocities already via Vx*Vx + Vy*Vy, even though we have not carried out the averaging 
transform mentioned before. The reason for this is that we used averaged classical guiding 
center coordinates (as given in the appendix) as zero order terms (with respect to h) in the 
general expansion ( p^ of the guiding center symbols. 

The term proportional to fi^ in (|77D specifies the lowest order quantum mechanical cor- 
rection to the classical guiding center Hamiltonian. It does not depend on the electrostatic 
field E = —'V(j). We will comment on its magnitude in the next section. 

Employing the same procedure to the more general case of a three-dimensional magnetic 
field and arbitrarily oriented electric field, the leading quantum correction to the classical 
Hamiltonian function turns out to be of second order in h as well. However, the results for 
the guiding center symbols and the guiding center Hamiltonian are notationally cumbersome 
and do not shed new light on our method. Therefore we refrain from writing them down 
explicitly in this paper. 
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4. SUMMARY AND CONCLUSION 



There are two major results of our investigations: First the product rule for gauge 
invariant Weyl symbols, equations (|36D , (|37|) , and (^Tj). They are a generalization of the 
well-known Moyal formula valid in the usual Weyl formalism. The leading order term in 
the ^-expansion p7|) of the Moyal bracket is proportional to the Poisson bracket, which is 
expressed in terms of derivatives with respect to position x and kinetic momentum k. The 
higher order terms in (^) are of a more complex structure and cannot be written as powers 
of the Poisson bracket operator. The question about their interpretation in terms of the 
modified phase space geometry in the presence of electromagnetic fields P2|-H3| may serve 
as an interesting starting point for further, more mathematical studies. 

The Weyl symbol Hw of the guiding center Hamiltonian, equation ([77D, represents the 
second major result of this paper. The method used to derive it makes extensive use of 
the product rule for gauge invariant Weyl symbols. The great advantage of this approach 
lies in the fact that the adiabatic parameter e can be incorporated into the gauge invariant 
Weyl calculus in a straightforward manner. Consequently, all expansions are carried out 
with respect to both e and Planck's constant h. 

Let us now investigate the importance of quantum mechanical effects on guiding center 
motion by comparing the magnitudes of the e^- and /i^-term in the guiding center Hamil- 
tonian ([77|) at low particle energies. Taking into account the quantization of energy levels 
we replace the gyration energies ^ BiV^ + Vy) with the corresponding harmonic oscillator 
eigenvalues {n + ^)hujB, where ujb = \qB\/{mc) is the cyclotron frequency at the position of 
the guiding center. In the absence of an electric field the guiding center Hamiltonian for a 
spinless particle takes the form (using vector notation) 

Hw = in+ l)hujB + ^'i"" (b^B -3\VB\^) + (bAB - \\/B\') , (79) 



with correct physical units restored in the way explained at the beginning of section 3. The 
second term in ( |75D is the adiabatic correction (with e set to unity) which now contains h due 
to energy quantization. Only for small oscillator (=gyration) quantum numbers n the lowest 
order classical and quantum mechanical corrections are of the same magnitude, otherwise 
the classical term dominates. From classical guiding center theory, however, it is known 
that at low particle energies the influence of adiabatic corrections on guiding center motion 
is negligibly small. Since the magnitude of the leading quantum mechanical correction in 
([f9|) does not depend on the particle energy, we conclude that quantum effects on guiding 
center motion may be neglected at all energy scales. Thus, when carrying out experiments 
with charged particles in inhomogeneous magnetic fields even at very low temperatures, the 
motion of the guiding center is described excellently by the lowest order classical equations. 
This is also true if an additional electrostatic field E is switched on, because according to 
([rzl) the magnitude of the lowest order quantum mechanical correction does not depend on 
E. As a general result of our investigations we may therefore say that guiding center motion 
is not effected by quantum mechanics. 



Note that if the quantum number n becomes too large in equation ([79|), the adiabatic 
correction dominates over the first (gyrative) term in the Hamiltonian. This parallels the 
breakdown of classical guiding center theory at large particle energies. 
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For reasons of completeness, let us finally compare our result (|77| ) for the guiding center 
Hamiltonian with the quantum mechanical calculation of Maraner, who investigated the 
motion of a charged spinning particle in a two-dimensional magnetic field. The interaction 
between the magnetic moment n of the particle and the external magnetic field is included 
into the Hamiltonian by the potential — /x • B. Since the component fiz of the magnetic 
moment parallel to B = B{x,y)ez is a constant of motion we may replace —fi B by the 
scalar term —fizB{x,y). Thus, the classical Hamiltonian function reads 

H = ^{vl + v^)-f^zBix,y). (80) 

The term —^zB{x,y) may be interpreted as a special case of a time-independent scalar 
potential y). A straightforward calculation shows that the corresponding guiding center 
symbol Hamiltonian is given by 

(fiAfi-3|Vfi|^) J( J-4/i,) -8^^5-2|y5|2- 

+ {BAB - I V5p) + 0{en\ e% e') , (81) 

where the magnetic field and its derivatives have to be taken at the guiding center position. 
Substituting fiz = —g^o's {g = gyromagnetic factor of the particle, (73 = ±1 for spin-| parti- 
cles), H{y takes the same form as the Hamiltonian operator derived by Maraner. However, 
in ( |8TD there are two expansion parameters e and h distinguishing between adiabatic and 
quantum mechanical corrections to the guiding center Hamiltonian whereas in [|l7l the only 
expansion parameter is the magnetic length Is = ^JhcJ{q\B\) . 
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APPENDIX: 

In this appendix we will briefly summarize the classical guiding center theory for the 
motion of a charged particle in a two-dimensional magnetic field B = B{x,y)ez and an 
electrostatic field E = —'V(f){x,y). We closely follow the Hamiltonian method developed by 
Littlejohn which employs non-canonical, but gauge invariant coordinates in phase space. 
For further details the interested reader is referred to |IO| , p]l|] . 

To incorporate the classical result into the symbol calculus we have to determine the 
Cartesian components Vx,Vy of the gyration velocity instead of the generalized gyrophase 
6 and magnetic momentum J. The reason lies in the difficulty of defining a quantum 
mechanical operator corresponding to the classical angle variable 6. This choice of guiding 
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center phase space coordinates has the disadvantage that in order to compute (X, Y, Vx, Vy) 



we cannot apply the elegant geometric method of [T^jni • 

Introducing the adiabatic parameter e in the standard way, i.e. by replacing the charge 
q by g/e, and scaling physical quantities according to equation (|53D we are looking for phase 
space functions (X, Y, 14, Vy) whose non-vanishing Poisson brackets are given by 

{v.^,v,} = \, {x.y} = ^. (Ai) 

Again, their different order with respect to e indicates the different time scales of motion. 
Making the perturbative ansatz 

oo 

Zi = Y,e''Z^[x,y,Vx,Vy) (A2) 

fc=0 

for the guiding center phase space coordinates {Zi\^^.^^ = (X, F, V^.,V^) we may calculate 
the coefficient functions Z^{x,y,Vx,Vy) order by order by solving the partial differential 
equations implied by the Poisson bracket relations ([Al|). To make sure that the functions 
(x, y, Vx, Vy) do not depend on the gauge, we have to write down the Poisson bracket in a 
gauge invariant form, 

{f^d} Qx dv dv dx ~^ e ^dv dvj ^ ^ '^'^ 

using only derivatives with respect to the physical (i.e. gauge independent) phase space 
coordinates of the particle, namely its position x and velocity v. 

The resulting guiding center coordinates are not uniquely determined by the Poisson 
bracket relations. In order to facilitate the following computations it is of advantage to choose 
them in the simplest form possible. To lowest order they are proportional to the related 
particle coordinates. Therefore, Vx,Vy are rapidly oscillating functions of time because the 
particle velocity components Vx,Vy depend on the gyration angle 9. 

Next we have to find the corresponding backward transformations, i.e. to write the 
particle coordinates {x,y,Vx,Vy) as functions of {X,Y,Vx,Vy). Inserting the result into the 
Hamiltonian function 



H 



the latter takes the form of an asymptotic series with respect to e with the coefficients being 
gauge invariant functions of the guiding center coordinates. To exclude rapidly oscillating 
terms in this expansion the Hamiltonian should depend on the gyration velocity components 
only by means of J = V^+Vy and its powers, because to lowest order in e J does not depend 
on the gyration angle 6. This can be achieved by carrying out a near-identity coordinate 
transformation {X,Y,Vx,Vy) — » {X ,Y ,Vx,Vy) to a new set of (averaged) guiding center 
coordinates. They have to satisfy the same Poisson bracket relations ( [Al| ) as the old ones 
so that the related magnetic moment J = V^ + Vy becomes a constant of motion and at the 
same time the dynamics of the guiding center position decouples from that of the gyration. 
Such a kind of symplectic transformation may be expressed in terms of Lie generators (cf. 

MM)- 
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As a result, up to second order in e the averaged guiding center coordinates read (leaving 
away the bar over them and using commas in the subscripts to denote differentiation with 
respect to the following coordinate (s)) 

^2 



X ^ X + €B~^Vy + 



{B^yV^-B^^Vy)Vy + O ( C'^ ) , 



Y ^ y - eB ^V:, - ^ (B^yV^ - B^^Vy)v^ + ^(e^), 
14 = B-^'-'v, + 



(A5) 
(A6) 



255/2 



By^ + B^^v^Vy + 2Byy + BB^y) 



+ 



16^9/2 



(-bBl - BB,,, + UBl - 5BB,yy) vl - 2 (4S,,S,, + 2BB,,y - ci) vlvy 

+ + + 23^2 - ^j,) + 2 {-UB^^B^y + 655^^, + ci) 

+ AB {B^xEx + BEx^x — 7B^yEy — SBEy^y ) Vx + 2B (c2 + 9B^xEy + 9B^yEx 



— 8BEx^y) Vy 



(A7) 



Vy = S-V2^^ _ __ (2By^ + B,yVxVy + Byy + BB,x) 



+ 



16^9/2 



2 {-UB^xB,y + GBB^xy - ci) + (235j - 1355 - 155f^ + 755,^,,) 
vlvy - 2 (AB,xB,y + 2BB,xy + ci) Vxvl + (l35j - 555,,, - 55f^ - 55,,,) 
- 25 (C2 - 75,5, - 75,5,) + 45 (-75 ,5, + 355,,, + 5,5, + 55,,,) t;, 

with the corresponding backward transformation 



X 



(A8) 



X 



y 



^ X - eB-'/'Vy 

= y + e5-^/V, 



252 



252 



(25,^2 + 5,,V;K; + B,xV; - 25,) + O(e^), 

(5,,v;^ + 5,y,y, + 25, if - 25,) + o(e^), 



(A9) 
(AlO) 



Vx = 5^2 1/ 



5 



(5,,y, + 5,,y,) V, - 5, 



+ 



1655/2 



(-115J + 55,,, - 35_2^ + 55,,,) l/J* - 4 (55,,5,, + 55,,, + ci/2) x 



VjVy + (5^4 + 555 - 155;^ - 355 ,,] VxVy' + 4 (5 ,5 , + 55 „ - ci/2) 



+ 4 {SB^xEx — BEx^x + B^yEy + BEy^y) Vx — 2{c2 — B^xEy — B^yEx) Vy 
Vy = 5V2 Vy + ^ [(5,,y, + 5,,y,) Vx - Ex] 



(All) 



+ 



1655/2 



4 (5,5, + 55,,, + ci/2) + (^_i5i^J _ ^bB^x. + B^ + 5BB,yy) x 
y,V, + 4 (55,5, + 55,,, - ci/2) VxVy' + (-W^ + 55,,, - 115,^^ + 55,,,) 
+ 2 (c2 + SB^xEy + SB^yEx — SBEx^y) V, + 4 [B^xEx + BEx^x + "^B^yEy — BEy^y) Vy 
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(A12) 



In ( |A7|) , (|A^), (|A11|) , (|A12|) , ci and C2 are arbitrary constants which remain uneffected by 
the requirement that the Hamiltonian has to be independent of the gyration angle up to 
second order in e (for more details concerning the ambiguity of guiding center coordinates 
cf. 10). In (|A9|)-( [XT2D the fields and their derivatives have to be evaluated at the guiding 
center position (X, F). 

In terms of the guiding center coordinates, the Hamiltonian reads 



H 



1 



16^2 L 



,y 



BB 



yy 



+ 4 (3E,i?,. - BEx,x + ^EyB,y - BEy,y) {v,' + Vy') - 8 {El + K 



(A13) 



which may be written in a more compact form by employing two-dimensional vector notation, 



1 e2 
H = -BJ + cb(X, Y) + 



BAB - 3\VBn + 4:{3E-VB - BV-E)J 



\E\ 
(AM) 



In (|A13|) and (|A14|) the fields and potentials have to be evaluated at the guiding center 
position {X,Y). For obvious reasons the last two expressions are denoted as the classical 
guiding center Hamiltonian. 
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